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APPLICATION TO HIGH-FREQUENCY COMBUSTION INSTABILITY 
by Charles E. Feiler and Marcus F. Heidmann 
Lewis Research Center 

SUMMARY 

The dynamic response characteristics of a lumped-element model of a gaseous - 
hydrogen injector were determined. The investigation was conducted to determine 
whether coupling between hydrogen flow oscillations and chamber pressure oscillations 
could be large enough to provide a mechanism for driving high-frequency combustion 
instability. The model was derived in terms of the injector design variables and was 
evaluated for three engines of different scales. The results are presented to show the 
effect of hydrogen density (temperature) on the coupling, since this parameter is useful 
in determining stability limits experimentally. A comparison of predicted and experi- 
mental stability limits was made which showed general agreement between experiment 
and theory. Such flow system coupling has been known to be of importance in low- 
frequency instability for some time. The present results show that, for gaseous flow 
systems, flow system coupling is an important factor in causing high-frequency instabil- 
ity also. 


INTRODUCTION 

An analysis was made of the dynamic response of the flow in a gaseous-hydrogen 
injector to determine whether the flow could couple with high-frequency chamber pres- 
sure oscillations to drive high-frequency combustion instability. While such coupling 
has been demonstrated for liquid flow at a low frequency, it apparently has not been con- 
sidered for gaseous flow at a high frequency. 

In this analysis, the hydrogen injector was modeled as a lumped-element system. 
Also, the equations describing the systems were linearized by a perturbation technique. 
Although the injector probably cannot be precisely described by a lumped-element 
system, the results should serve to illustrate the importance of flow system coupling on 
instability. 

The analysis is presented in terms of the injector and engine parameters that are 



ordinarily varied in engineering design. It has been examined over a range of parame- 
ters by considering three engines of widely differing scales. The results are presented 
to show the effect of increasing hydrogen density (decreasing hydrogen temperature) on 
the coupling between oscillations in hydrogen flow and chamber pressure. Hydrogen 
density was chosen for the independent parameter to facilitate comparison of the results 
with the experimental data of reference 1. These experimental studies (ref. 1) employed 
a hydrogen temperature -ramping technique in which the hydrogen injection temperature 
was gradually decreased until instability occurred. This transition temperature was 
then used to rate the various injector schemes and engine designs; the lower the transi- 
tion temperature, the more stable the system was considered to be. The comparison of 
these data with the present analytical results gives some indication of the significance 
of the analysis in a qualitative manner. 


SYMBOLS 

2 2 

A cross-sectional area, in. (cm ) 

C capacitance term (eq. 19(a)), sec 

g gravitational constant, 32. 174 lb mass-ftyib force-sec 2 (1 kg - m/N - sec 2 ) 

I inertia term (eq. 19(a)), sec 

L hydrogen annulus length, in. (cm) 

(L/A) eff effective value of (L/Aj), in. "* (cm - *) 

M hydrogen mass in dome, lb mass (kg) 

N response factor (eq. (1) or (2)) 

O/F oxidant -fuel mixture ratio 

9 9 

P c chamber pressure, lb force/in. (N/m ) 

2 2 

P^ dome pressure, lb force/in. (N/m ) 

2 2 

Pj annulus entrance pressure, lb force/in. (N/m ) 

9 9 

Pg orifice entrance pressure, lb force/in. (N/m ) 

Rp Rg pressure -ratio functions, (eq. (22)) 
s Laplace operator 

T temperature, °R (°K) 

t time, sec 

U velocity, ft/sec (m/sec) 
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V dome volume, in. ** (cm**) 

W flow rate through injector, lb mass/sec (kg/sec) 

Wjj flow rate of burned hydrogen gases, lb mass/sec (kg/sec) 
W^ n dome inlet flow rate, lb mass/sec (kg/sec) 

Wj. total flow rate, hydrogen and oxygen, lb mass/sec (kg/sec) 

y specific -heat ratio 

9 phase angle between P^, and W* 

p hydrogen density, lb mass/ft (kg/m ) 

time delay constant for burning, sec 
p reactance, eq. (26) 

o> frequency, rad/sec 

Subscripts: 
lox liquid oxygen 

max peak amplitude of sine wave 
noz nozzle 

1 hydrogen annulus 

2 hydrogen orifice 
Superscripts: 

average value 

' perturbation quantity, (x - x)/x 


ANALYSIS 

Conceptual Philosophy 

The analysis was conducted to determine the extent to which the hydrogen flow could 
couple with chamber pressure oscillations to drive them. The analysis differs from the 
usual frequency-response analysis in that the frequency and oscillation mode character- 
istics were fixed by engine geometry, and the effects of other parameters on the coupling 
were examined. The result of the analysis is a response factor N that gives the real 
part of the hydrogen-flow-rate oscillation in phase with the chamber pressure oscilla- 
tion. The response factor as derived in reference 2 for sinusoidal oscillations is 
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( 1 ) 


„ „„„ a 

N = cos 0 

(P r ) 

' c, max 7 

where P’ = (P* „, 0 „)sin wt and W* = W* sin(o>t + 9 ). 
c c j max max 

Before the response factor for the hydrogen flow is calculated, a brief discussion of 
its significance and use is presented. A system with a positive response factor (in-phase 
coupling) should provide a driving source, and a system with a negative response factor 
(out-of -phase coupling) should provide damping. In this isolated system, neutral stabil- 
ity should correspond to a response factor equal to zero. In a rocket engine, stability is 
controlled by several processes that act simultaneously and may be sources of either 
driving or damping. In the real case of interest then, stability would occur if each 
response factor were zero; however, a more realistic criterion would seem to require 
only that the sum of the response factors be zero or that the gains of the entire system be 
equal to the losses. This criterion holds for a volume in which the pressure oscillates 
uniformally. It is also applicable to a transverse mode in a rocket combustor with uni- 
form injection and a uniformally distributed nozzle. The following expression defining 
stability in terms of the response factor and representing the foregoing criterion was 
assumed: 


= — N + ^N + ^N =0 

W t noZ W t lox W t H 2 


( 2 ) 


where the response factors N noz , Nj qx , and N H are for the nozzle flow, the oxygen 

2 

vaporization or combustion, and hydrogen flow, respectively. In the summation, each 
response factor was weighted by the fraction of the mean total weight flow involved in the 
particular process. To use this relation, an evaluation of the oxygen and nozzle response 
factors is needed. For oxygen, the peak linear response factor Nj ox of 0. 55 from 
reference 2 was used as an approximation. 

For the nozzle response factor, the simplest assumption was to consider the nozzle 
as a purely resistive element and to neglect velocity effects. From the continuity rela- 
tion, 


W = -PAU 


which, in perturbation form, gives 



( 3 ) 


( 4 ) 
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where P = P y and P ’ = yp'. 

According to equation (1) the nozzle response factor for this model is 

N „ = - 1= -0.833 for y = 1.2 
noz y 

An approximation to the value of the hydrogen response factor at the stability limit 
becomes possible through the use of equation (2) and these values of the nozzle and oxygen 
response factors. 


Hydrogen System Model 

In figure 1, a single coaxial injector element of the hydrogen flow system is shown. 
The coaxial injector type is the predominant choice in practice for the hydrogen-oxygen 
propellant combination. It consists of an annular orifice (area Ag) connected by an 
annular passage (area Aj and length L) to a supply dome or cavity (volume V). It is 
assumed that the dome is supplied at a constant flow rate W^ n - A transfer function for 
these elements that dynamically relates the flow rate to the chamber pressure is first 
derived; from this the response factor can be derived directly. As stated earlier, the 
equations in perturbation form are linearized. 

Supply dome (capacitance) . - For the fluid in the dome, a mass balance gives 

— = W in - W (! 

m m 


Since W^ n is a constant and M = pV, in perturbation form equation (5) becomes 


PgV 

w 


sp T = - W' 


( 6 ) 


where s is the Laplace operator. Throughout this development, the adiabatic relation 
between pressure and density is used. In perturbation form this relation is 

p’ = yp’ where y = y(T) (7) 

Substitution of equation (7) into (6) yields 
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( 8 ) 


El 

y W 



Annular entrance loss (resistance). - The flow from the dome into the annulus is 
represented by a resistive entrance loss for which the orifice flow relation is assumed: 

W = KA ± V(P d - P^Pg (9) 

With the aid of equation (7), in perturbation form, 

where AP^ = P^ - P^ (10) 

Annulus flow loss (inductance) . - The flow in the annulus is assumed to be of slug 
form. The associated inertia is obtained from the momentum and continuity relations 
with the assumption of an average density and velocity over the length of the annulus. 


w' = I 


± p’ + (I -ZLjPj 
APi d \y AP t ' 1 


P 1 - P 2 = 


L dW 
gA 1 dt 


or 


P l P i - P 2 P 2 * ~ sW ' 

gAl 


If it is assumed that Pj = P 2 , then 


P’ - P, = -^ L sW' 

g P 2 A l 


Final orifice (resistance) . The final orifice yields 

W = KA 2 ys* 2 - P c )Pg 


( 11 ) 


( 12 ) 


(13) 


( 14 ) 


or 
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where APg = P 2 " P c 


(15) 


W = 


1 

2 




Burning process . - In addition to the relations presented thus far, it is assumed that 
before the hydrogen can contribute to the chamber pressure it must, in some way, by 
diffusion and mixing and reattion, be brought to the chamber temperature. This step is 
represented simply by a delay time. Thus, 

W b (t) = W(t - r b ) (16) 


and 


TkS 


W b = W'e b 


(17) 


The perturbation expressions (eqs. (8), (10), (13), (15), and (17)) may be combined 
to give the transfer function 



where 


and 


C 


= £V 

y W 



~\ 


> 


J 


(19) 


Following the procedure of reference 2, the response factor N given by equation (1) is 
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where 



and 



By defining two dimensionless parameters, 




and 


*2 S 



( 21 ) 


( 22 ) 


( 23 ) 


equations (21) and (22) may be more simply expressed: 



1 


( 24 ) 


W 

b, max 



Co? 

A final simplification is made by defining a reactance i f/ where 



This results in the following modified response factor: 


N 



+ 




(25) 


(26) 


(27) 


ANALYTICAL RESULTS 
General Behavior of Model 

Equation (27) represents analytically the dynamic response of the flow for the injec- 
tor element of figure 1. This equation is plotted in figure 2 for both positive and negative 
values of if/ and for values of from 0 to n radians. In this figure, the response 

of the flow system alone is given by the curve for cor^ = 0. For = 0, the hydrogen 
flow system tends to be a damping process, as shown by the negative response factor, 
because the flow oscillations are out of phase with chamber pressure oscillations. At 
large if/, either positive or negative, the response factor approaches zero because the 
flow no longer responds to pressure oscillations and because the phase angle approaches 
90°. At the unique point if/ - 0, the response curve attains its maximum value (negative). 
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At this point, the system is, in a sense, tuned and the flow oscillations are exactly 180° 
out of phase with pressure oscillations. The value of the response factor then depends 
only on the resistive pressure drop terms. 

At = 0, W* and are 180° out of phase as just stated. As is increased 

from 0 to ir/2, the response factor becomes less negative and is zero at cor^ = ir/2, 
where W' and P ' are phase displaced by 90°. As is increased from -n/2 to it, 

the response factor becomes more positive and attains a maximum at cor^ = u. At this 
point, W* and P^ are exactly in phase. The further increase of from v to 2n 

is not shown in figure 2, but it results in the response factor decreasing to its value at 
u)T b = 0. A similar behavior occurs at other values of ip, but it is more obscure be- 
cause of the effect of xp on both w\ _.__./P' and 6. The limiting curve of the 

response factor is given by the dashed line. This curve is obtained at each value of xp 
by adjusting cor^ so that 9 = 0. 


Parametric Study 

The response factor was calculated for three engines to show the effect of the dif- 
ferent parameters. These were the engine of reference 1, which is similar to the RL-10 
and engines A and B, about the size of the M-l and J-2, respectively. The three engines 
differ greatly in scale, as shown by the parameters listed in table I. Each of the param- 
eters Ag, V, and L/A was evaluated at one-half and twice the value listed in table I 
to show the sensitivity of the response factor to the parameters. 

For the inertia term, an effective value of (L/Aj), (L/A) e j f composed of the sum of 
two contributions, was used. The contributions considered were those from the final 
orifice and the annular passage. No correction for effective length was made in this 
evaluation. The individual injector elements were assumed to behave as isolated units. 
The volume for an individual element was approximated as the ratio of the volume of the 
dome to the number of elements. Since the calculations were made with the use of total 
hydrogen flow rates and injection areas, the total dome volume could be used in the 
numerical procedure. 

The time delay was unknown but its use is justified as a simple representation 
of the history of the hydrogen gas between injection and combustion. For a typical hydro- 
gen injection velocity of 500 feet per second (152 m/sec) and a of 0.0001 second, the 
events represented would have to occur in a space of 0. 6 inch (1. 52 cm). This distance 
approximates the standoff distance that has been observed for the flame; however, the 
significance of may be more complex or even unrelated to such a simple concept. 

The value of (0. 00009 sec) used for the engine of reference 1 more closely repre- 
sented the data, as is shown later. 
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TABLE I. - CALCULATION PARAMETERS 


Parameter 


Engine 



Reference 1 

B 

A 

Chamber pressure, P , 

300 (20. 7xl0 5 ) 

680 (46.9X10" 5 ) 

1000 

lb force/in. 2 (N/m 2 ) 



(69.0X10 5 ) 

Total flow rate, Wrp, 

65 (29.5) 

463 (210) 

3420 (1550) 

lb mass/sec (kg/sec) 




Frequency, o>, rad/sec 

21 350 

11 840 

5340 

Oxidant -fuel ratio, O/F 

5.5 

5.5 

5.5 

Hydrogen orifice area, 

4.62 (29.8) 

25. 7 (166) 

69 (445) 

Ag, in. 2 (cm 2 ) 




Dome volume, V, 
in. 3 (cm 3 ) 

46.0 (754) 

274 (4490) 

1690 (27 700) 

Hydrogen annulus length 

0.0677 (0.0267) 

0.0102 (0.0040) 

0. 0046 (0. 0018) 

to area ratio, (L/A)eff, 
in. (cm -1 ) 




Hydrogen annulus area, 

9. 12 (58. 8) 

68. 2 (440) 

205.5 (1325) 

Aj, in. 2 (cm 2 ) 




Time delay constant for 
burning, sec 

0. 00009 

0.0001 

0.0001 

Hydrogen density, p, 

0. 1 to 4. 04 

0. 1 to 4.04 

0.1 to 4.04 

lb/ft 3 (kg/m 3 ) 

(1.60 to 6.47) 

(1.60 to 16.47) 

(1.60 to 6.47) 

Temperature, T, 

a ~645 to ~50 

b ~1400 to ~55 

c ~2150 to ~60 

°R (°K) 

(358 to 27. 8) 

(778 to 30. 6) 

(1195 to 33. 3) 


a At 20 atm (20. 3xl0 5 N/m 2 ). 
b At 43 atm (43.6xl0 5 N/m 2 ). 
c At 70 atm (71xl0 5 N/m 2 ). 


Since the Mach number of the flow was low, one further approximation that was made 
was the use of a constant average density in calculating the mean values of the pressure 
drop terms APj and APg. A discharge coefficient of 1, was also assumed. 

In figure 3, the variation of the response factor with hydrogen injection density is 
shown with the injector design parameters as variables. The effect of ramping the hydro- 
gen temperature is readily seen. Figure 3(a) shows the oxidant -fuel ratio O/F as a 
variable; increasing O/F tends to increase the magnitude of the response factor. Fig- 
ure 3(b) shows a similar result with the final orifice area A 2 as a variable. The major 
effect is upon the resistive terms. Increasing either O/F or A 2 at a constant total 
weight flow W t causes a reduction in the pressure drop AP 2 across the orifice A 2 . 

As a result, the coupling between W* and P' is greater as is the absolute value of the 

V 

response factor. For the A and B engines, this coupling was out of phase while for 
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II 


the reference 1 engine it was in phase. 

Figures 3(c) and (d) show the effect of V and (L/A)^, respectively. The effect of 
increasing either term is qualitatively the same, since an increase in either results in a 
change in the value of \f/ to more negative values. Neither V nor (L/A) e ^ affects the 
resistive terms. The result is that the peak or any point on the i// curve (fig. 2) occurs 
at a different density and at a different resistance, which shifts the response factor 
curves so that similar points, such as the minimums and zero response factor, occur 
at different densities. In figures 3(a) and (b), where the primary effect was on the 
resistive terms, these corresponding points occur at about the same density for a partic- 
ular engine. 

The general effect of increasing density, as shown by all the curves in figure 3, was 
to increase the absolute value of the response factor. Increasing the density has two 
primary effects: first is that of lowering the resistance, or pressure drop, which re- 
sults in a greater coupling, and second is that of shifting the xj/ toward more negative 
values. The major part of the effect on ^ occurs by the increase in the capacitive 
term C, which is directly proportional to the density. The inertia term I is not strongly 
influenced by density. 

One other general point can be made about the nature of the curves of figure 3. The 
density at which the response factor is zero occurs when the total angle 6 is n/2 . This 
density differs among engines because the resonant frequency and the i// values differ. 

The parameters that were evaluated are related to injector design and to the scale 
of the engine. Other parameters such as contraction ratio would also influence the 
response factor by affecting the chamber pressure, weight -flow rate, or resonant fre- 
quency. These parameters have not been included as a part of this study. 

The results of the analysis show an appreciable coupling of the hydrogen flow oscilla- 
tions with the chamber pressure oscillations. Hydrogen density is shown to be an im- 
portant factor in the coupling. For a given engine, the variation of selected injector 
parameters over a range of 4 to 1 did not change the basic character of the response 
factor for that engine. This basic character appears to be determined by the scale of the 
engine, in that one of the important parameters was the resonant frequency. For a con- 
stant time delay, as was assumed, the phase angle 9 is strongly influenced by the con- 
tribution of to 6. This influence, in turn, determines whether the basically out - 
of -phase coupling between flow and pressure oscillations can be shifted to become in 
phase. 


COMPARISON OF EXPERIMENT AND ANALYSIS 

The engine of reference 1 offers an opportunity for comparing experimental and 
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analytical results. The hydrogen-temperature stability limit of this engine was experi- 
mentally determined as a function of O/F for three hydrogen injection areas with the 
same oxygen injector. Figure 4 shows the response -factor - density curves calculated 
for this engine over a range of O/F values and for three orifice areas Ag. The 
response factor increased with density for this particular engine. As a result, decreas- 
ing hydrogen temperature (increasing density) should cause the engine to become ion- 
stable. The value of the hydrogen response factor at the stability limit was calculated 
by using equation (2) and the values of nozzle and oxygen -system response factors given 
previously. From figure 4, the density corresponding to the hydrogen response factor 
may be obtained. 

In figure 5, the stability limit data of reference 1 are replotted to show hydrogen 
density rather than temperature as a function of O/F. The data show the hydrogen den- 
sity separating stable and unstable operation. The curves are predicted from the anal- 
ysis. The curve A 2 = 4. 62 square inches (29. 8 sq cm) was obtained by selecting a value 
of the time delay that reasonably approximated the data. For this purpose, was 
0. 00009 second, and this value was used for the remaining injectors. With this proce- 
dure, the analysis overpredicted the injector of Ag = 1. 41 square inches (9.1 sq cm) by 
about 50 percent. The injector with A 2 = 0. 77 square inch (5. 0 sq cm) was essentially 
stable at the lowest density attainable in the experiment. The calculation for this in- 
jector also indicated that it was stable to a density greater than that attainable. 

From this comparison with the data of reference 1, it can be seen that, while agree- 
ment was not perfect, the trends with both O/F and the hydrogen injection area were 
predicted. Better agreement could probably have been obtained by allowing the time 
delay to vary with the hydrogen injection area; however, too little is known about the 
time delay to justify its variation. In making the comparison, it was also assumed that 
changing the hydrogen injection area did not influence the oxygen response factor. Relax- 
ing this assumption might also improve the agreement. The important conclusion, how- 
ever, can still be made: the hydrogen flow is sufficiently sensitive that the coupling be- 
tween hydrogen flow oscillations and chamber pressure oscillations can drive instability. 

Similar experimental data were not available for other engines. The present anal- 
ysis indicates that engine A (M-l size) should be a stable configuration with regard to 
the hydrogen flow system, as shown in figure 3(a), where the response factor is gen- 
erally negative. For engine B (J-2 size), figure 3(a) shows that the response factor in- 
creases rapidly at a density greater than 2. 5 pounds per cubic foot (40 kg/cu m). At the 
pressure level of the B engine, this density occurs at about 60° R (33° K). As more 
data become available, they may provide insight into the importance of various param- 
eters in the model that has been presented and thus permit a more realistic model to be 
derived. This method of analysis also predicts qualitatively the increased stability (lower 
hydrogen temperature or higher density) of large oxygen jets reported in reference 1. 
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In this case, a lower value of oxygen response factor for large jets and drops as predicted 
in reference 2 is used in the summation of response factors (eq. (2)). Such a trade off 
between oxygen and hydrogen response factors gives higher density (lower temperature) 
stability limits for larger diameter oxygen jets. 


SUMMARY OF RESULTS 

A lumped-element analysis of the response of a hydrogen-feed system to high- 
frequency chamber -pressure oscillations was made. The results may be summarized 
as follows: 

1. For the injector model assumed, there was sufficient coupling between flow and 
pressure oscillations for this coupling to be an important factor in instability in gaseous - 
hydrogen systems. 

2. Hydrogen density, an important variable in experimental studies, was also shown 
to be important in the analysis. Coupling of the flow and pressure oscillations increased 
with density. 

3. Although appreciable coupling of the flow system to pressure oscillations was 
found, it was basically an out-of -phase coupling. The representation of the combustion 
process by a simple delay time concept changed the coupling from out of phase to in 
phase or from damping to driving. 

4. A general agreement in trend between available experimental data and analytical 
predictions was found. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 11, 1967, 

128-31-06-02-22. 
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Figure L - Hydrogen injector element. 
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Figure 2. - Generalized response factor. 
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(a) Effect of oxidant-fuel ratio. 


Figure 3. - Variation of 
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(b) Effect of hydrogen orifice area. 
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(c) Effect of dome volume. 



2 4 6 8 10 20 40 60 

Hydrogen density, p, kg/m 3 


645 322 145 104 87 78 66 50 

Hydrogen temperature (at 300 psi), T, °R 


358 179 80.6 57.8 48.3 43.3 36.7 27.8 

Hydrogen temperature (at 20. 7xl0 5 N/sq m), T, °K 

(d) Effect of inertia. 


Figure 3. - Concluded. 
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Figure 4. - Variation of response factor with hydrogen density for engine of 
reference 1. 
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Hydrogen density, p, lb mass/ft 



Figure 5. - Comparison of experimental and calculated hydrogen-density stability 
boundaries. 
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